Radiative capture reactions in lattice effective field theory 
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We outline a general method for computing nuclear capture reactions on the lattice. The method 
consists of two major parts. In this study we detail the second part which consists of calculating 
an effective two-body capture reaction on the lattice at finite volume. We solve this problem by 
calculating the two-point Green's function using an infrared regulator and the capture amplitude to 
a two-body bound state. We demonstrate the details of this method by calculating on the lattice the 
leading Ml contribution to the radiative neutron capture on proton at low energies using pionless 
effective field theory. We find good agreement with exact continuum results. 
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Introduction. — A long standing goal of nuclear physics 
and astrophysics is to understand the nuclear reactions 
that made the elements of nature, from the lightest nu- 
clei in the early universe to the light, medium, and heavy 
nuclei synthesized in stars. In the past few years there 
has been progress in ab initio calculations of nuclear reac- 
tions which goes beyond transitions from one bound state 
to another. Recent work includes calculations using the 
No-Core Shell Model and resonating group method PQH], 
Fermionic Molecular Dynamics [3|, the coupled-cluster 
approach [4], and variational Monte Carlo 0. There 
has also been some progress using lattice calculations in 
finite periodic volumes to analyze coupled-channel scat- 
tering [5H5] and three-body systems [lOj. However there 
is no general formalism for calculating inclusive and ex- 
clusive reactions from lattice simulations. In this letter 
we present the first steps towards a general method for 
computing nuclear reactions on the lattice. Our discus- 
sion here focuses on radiative capture reactions using the 
formalism of lattice effective field theory. Our results will 
also have direct applications to ultracold atomic systems 
where there are interesting phenomena analogous to nu- 
clear capture reactions. In this case radio-frequency pho- 
tons are used to induce the association of weakly-bound 
molecules [TItiTS] . 

Lattice effective field theory combines the theoretical 
framework of effective field theory (EFT) with numerical 
lattice methods. A review of lattice effective field theory 
calculations can be found in Ref. [IS]. The method has 
been applied to nuclei in pionless EFT [lT| and chiral 
EFT [IMSQ]. In Ref. [21] a new technique was developed 
which allowed for general wavefunctions to be used as 
initial and final states in the lattice EFT calculations. It 
has since been realized that cluster wavefunctions can be 
used to study the scattering and reactions of continuum 
states. The general strategy involves separating the cal- 
culation into two parts. The first part of the method is to 
use projection Monte Carlo to determine a multi-channel 



adiabatic lattice Hamiltonian for the participating nuclei. 
For example, let us consider a two-body capture reaction 
involving three nuclei, 1-f 2 — > 3 + 7. Let |ri2) be a initial 
cluster wavefunction with separation vector ri2 between 
nuclei 1 and 2. We consider all possible values for ri2 and 
use Euclidean time propagation to construct projected 
states, |ri2)j = exp (— |ri2). For large t, the set of 
states |ri2)( will approximately span the linear space of 
continuum states for nuclei 1 and 2. Using these states 
we can calculate one-photon transition matrix elements 
with nucleus 3. The details of this projection technique 
and construction of the multi-channel adiabatic lattice 
Hamiltonian will be discussed in forthcoming publica- 
tions. The focus of this letter is describing how to do 
the second part of the reaction calculation. This second 
part uses the adiabatic lattice Hamiltonian for the par- 
ticipating nuclei to calculate nuclear reaction rates. For 
the case of nucleon-nucleus and nucleus-nucleus capture 
reactions, the problem is equivalent to the radiative cap- 
ture of two point particles forming a two-body bound 
state. In this letter we show how this calculation is done 
by using the well-known example of radiative capture of 
a proton and neutron to form a deuteron. In our analy- 
sis we use pionless EFT at leading order. We will show 
that the continuum results which have been obtained in 
Ref. [52, 23 can be accurately reproduced on the lattice. 

Neutron capture. — The radiative capture process 
p{n,j)d is important in nuclear physics as it provides 
stringent bounds on the primordial deuterium abun- 
dance. This abundance is sensitive to the amount of 
baryonic matter in the universe. At very low energies, 
this reaction proceeds primarily through the Ml transi- 
tion from the incoming p + n spin-singlet s-wave state 
to the final spin-triplet s-wave deuteron bound state. It 
is known that these initial and final s-wave states are 
non-perturbative. Therefore calculating the Ml transi- 
tion using lattice EFT provides a non-trivial check of the 
method developed here. We add that at the energies rel- 
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evant in Big Bang nucleosynthesis, p{n, is dominated 
by both the Ml and El transitions. However, the initial 
state p-wave interaction associated with the El transi- 
tion is perturbative. So while we do not include the El 
contribution in the lattice EFT calculation described in 
the following, this contribution can be added trivially. 

The dominant Ml contribution can be calculated with 
the leading-order Lagrangian, 
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where D^^ = 9^ + ie(l -I- T3)A^/2 is the covariant deriva- 
tive, Bi the magnetic field and Ki — 2.35 the isovector 
nucleon magnetic moment. We take the nucleon mass as 
M = 939 MeV. The Pauli matrices ai act on the spin in- 
dices and Ti act on the isospin indices of the nucleon field 
N. The Pauli matrices a2T2Ti and a2(JiT2 project the in- 
coming and final state nucleons onto the spin-singlet and 
spin-triplet channels, respectively. 



iAa{p) 



FIG. 1. Elastic scattering of the nucleon. The index a — s,t 
corresponds to the singlet and triplet channels. The "• • • " 
represents iteration of the interaction. 

The elastic scattering diagrams in Fig. [l] give the am- 
plitude 
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where A is the renormalization scale and the index a = 
s,t corresponds to the spin singlet and triplet channels, 
respectively. In the spin-singlet channel, we tune the 
EFT couplings in the continuum such that there is a 
shallow virtual state with a scattering length a = —23.7 
fm with 4:Tr/{Mcs) -l- A = 1/a. In the spin-triplet chan- 
nel we tune the EFT couplings to form a shallow bound 
state with binding momentum 7 = 45.7 MeV where 
An/iMct) -f A = 7. 

The leading order contribution to the capture ampli- 
tude comes from the diagrams (a) and (6) in Fig. [2] Using 
the center-of-mass (cm.) kinematics with p the neutron 
momentum and k the photon momentum, we get |22l 123] 





FIG. 2. Radiative capture diagrams for p{n, 7)d . Wavy line 
represents the photon, dashed line the deuteron, and blob the 
set of possible initial state interactions. 



in the continuum 
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Un represents the spinor wave functions for the nucleon 
fields, and Z is the wave function renormalization factor 
defined as the residue at the deuteron propagator pole. 
In the lattice calculation, we will reproduce the reduced 
matrix amplitude Aici a-nd not concern ourselves with 
the overall factors associated with the standard normal- 
ization of non-relativistic cross sections. 

The amplitude A^c is related to the expectation value 
('i/'B |Oem of f relevant electromagnetic operator be- 
tween final bound state wave function 4'b{i') and initial 
state wave function ipi{r). On the lattice, it is straightfor- 
ward to calculate the bound state wave function 
from the discretized Hamiltonian. The incoming wave 
function ipi{r) is problematic since we work with peri- 
odic finite volumes. One could do a finite volume analy- 
sis as in Ref. [8 , but this mixes in the elastic phase shift 
with the capture amplitude. Here instead we calculate 
the inelastic process on the lattice using the interacting 
two-point retarded Green's function. We write, in lattice 
units, 



G{E;x,y) = {y\ 
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where G{E] x, y) is the Green's function for propagation 
from X to y with the strong interaction Hamiltonian Hg 
in the incoming spin-singlet channel. The bound state 
wave function is normalized such that = 
l/(87r7). For on-shcU incoming particles, the pre- factor 
— E — ie is just —ie in Eq. (jij), as prescribed by the 
Lehmann-Symanzik-Zimmermann reduction formula for 
the scattering matrix. 

For the lattice calculation of Ai in Eq. Q we use the 
Hamiltonian lattice formalism with continuous time. The 
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lattice Hamiltonian is derived from discretizing the con- 
tinuum Hamiltonian corresponding with Eq. ([T]) . We use 
the simplest possible discretization with nearest neighbor 
hopping terms for the kinetic energy term and single- 
site two-body contact interactions. We use a lattice 
spacing oi b = 1/100 MeV^^ for scattering momenta 
p < 40 MeV and b = 1/200 MeV"^ for scattering mo- 
menta p > 40 MeV. In each case we use periodic bound- 
aries in all spatial directions and consider cubic lengths 
up to L = 2406. The lattice spacing b and box sizes L 
are chosen so that 27r/L ^ p <C 7r/& for each momentum 
P- 

The strength of the triplet-channel coupling q on the 
lattice is tuned to generate a bound deuteron with en- 
ergy —2.225 MeV. This sparse matrix lattice eigenvec- 
tor calculation also gives the bound state wave function 
tpBif)- For the spin-singlet channel we use Liischer's 
finite volume method |24| to determine the interaction 
coefficient c^. Liischer's finite volume formula relates the 
two-particle energy levels in a periodic cube to the scat- 
tering phase shift. We use Liischer's method to tune the 
lattice coupling to produce the neutron-protron scatter- 
ing length, a = —23.71 fm. 

In the continuum limit we can analytically calculate 
the dependence on the parameter ie at infinite volume. 
The exact result is a generalization of Eq. ^ 
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with pe = \/ p^ -\- iMe. The parameter e serves as an in- 
frared regulator which exponentially suppresses the con- 
tribution from scattering at large distances. This is a key 
fact which we use to remove finite volume errors from the 
lattice calculations of A4 in Eq. For nonzero e the 
finite volume error becomes exponentially small as a func- 
tion of box length L. This allows us to perform lattice 
calculations for values of L where the finite volume error 
is negligibly small. We then use the fact that the scat- 
tering amplitude is smooth in e to extrapolate to e = 
using a simple linear extrapolation. We note that there 
are computational methods such as the Lorentz Integral 
Transform method |25| which also compute the ampli- 
tude as a function of e. However the simple extrapola- 
tion method appears to be the most efficient approach 
for the two-body capture problem we analyze here. A 
similar complex energy method was used in Ref. [26] for 
Faddeev-Yakubovski calculations. 

For convenience we define the parameter S = eMp^ . In 
our lattice calculations we compute the scattering ampli- 
tude at volumes where the finite volume error for our 
chosen value of e is negligible. We then extrapolate the 
lattice data to the ie — ?> «0+ limit using the linear fit 
So + si S. This extrapolation corresponds with taking the 
infinite volume limit. The extrapolation for several val- 
ues of p are shown in Fig. [Sj In the top panel we show 



results for \ Ai\ which is the ratio of to the contin- 

uum limit result at 6 = 0, |A^c(0)|. In the bottom panel 
we show results for which is the ratio of arg[Ai{S)] to 
the continuum hmit result at S — 0, arg[Mc{0)]- 

The errors of the extrapolation fit are indicated by the 
So values in the legends in Fig. [3[ The deviation of the 
extrapolated value of |A^| and from the value 1 is a 
measure of the lattice spacing error. The extrapolated 
lattice numbers agree with the continuum results for \M\ 
to within 2% and for to within 4%. We don't present 
results for other lattice spacings, but the comparison with 
the exact continuum limit results makes it clear that the 
lattice discretization errors are under good control. 
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FIG. 3. Linear extrapolation of lattice results (solid curves). 
The top panel shows_|A?| = \M{S)\/\Mc{5 = 0)|, and the 
bottom panel shows <jf> — arg[M(<5)]/ arg[Afc((5 ~ 0)]. 

In Fig. |4] we show the lattice results from Eq. 0) in 
comparison with the continuum results from Eq. (IsF for 
5 = 0.6, 0.4 and results extrapolated to (5 = 0. In the top 
panel we show \A4\ and in the bottom panel we show the 
phase angle in degrees. We see that the lattice results 
reproduce the continuum results in all cases with errors 
no more than a few percent. This residual error can be 
attributed to lattice spacing discretization effects. 

Discussion. — We have proposed a general formalism 
for ab initio calculations of radiative capture cross sec- 
tions on the lattice. The first part of method requires 
performing lattice simulations to determine the matrix el- 
ements of a multi-channel adiabatic Hamiltonian, reduc- 
ing the incoming state to an effective two-body system 
and the outgoing nucleus to a single body. The second 
part of the method uses the two-point Green's function 
for this adiabatic Hamiltonian to calculate the radiative 
capture cross section. 

The numerical details of the first part of the method 
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FIG. 4. Comparison of lattice results with continuum results 
from Eq. ([5| for 5 = 0.6, 0.4 and results extrapolated to 5 = 0. 
In the top panel we show \A4\ and in the bottom panel we 
show the phase angle (fi in degrees. 



will be discussed in forthcoming publications. In this 
study we have demonstrated the second part of the 
method using the example of radiative neutron capture 
reaction p{n,j)d at leading order in pionless EFT. We 
have shown that extrapolations in the parameter S are 
an effective way to remove finite volume errors in the lat- 
tice calculations. We were able to reproduce the known 
continuum result for the capture amplitude at infinite 
volume with an error of within 2% for the magnitude 
and an error of within 4% for the phase angle. These 
remaining errors are consistent with the size of the small 
discretization errors due to lattice spacing. 

This method should be applicable for calculating pho- 
tonuclear reaction rates involving many nuclear systems. 
This includes halo nuclei with a single valence nucleon. 
For example, in ^^C(7i, 7)^'^C the carbon-15 nucleus is 
a halo nucleus with a neutron separation energy of 
only 1.22 MeV. At low energy the incoming state can 
be treated as a point-like neutron and a tightly-bound 
carbon- 14 core in halo EFT jJT]. But ab initio lattice 
calculations should be able to go well beyond halo sys- 
tems. This method could also be applied to more com- 
plicated systems such as alpha capture on carbon-12 or 
oxygen- 16 as occurs in helium burning. We are hopeful 
that such systems will studied in the near future using 
methods outlined here. 
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